Nondestructive volumetric imaging of tissue 
microstructure with benchtop x-ray phase- 
contrast tomography and critical point drying 

Adam M. Zysk, 1 Alfred B. Garson III, 2 Qiaofeng Xu, 2 Eric M. Brey, 1 Wei Zhou, 1 
Jovan G. Brankov, 3 Miles N. Wernick, 1 ' 3 ' Jerome R. Kuszak, 5 ' 6 and 
Mark A. Anastasio 2 '* 

'Department of Biomedical Engineering, Medical Imaging Research Center, Illinois Institute of Technology, 

3440 S. Dearborn St., Suite 100, Chicago, Illinois 60616, USA 
"Department of Biomedical Engineering, Washington University in St. Louis, One Brookings Drive, St. Louis, 

Missouri 63130, USA 

'Department of Electrical and Computer Engineering, Medical Imaging Research Center, 
Illinois Institute of Technology, 3440 S. Dearborn St., Suite 100, Chicago, Illinois 60616, USA 
'Department of Ophthalmology and Visual Sciences, University of Illinois at Chicago, 1855 West Taylor Street, 

Chicago, Illinois 60612, USA 
'Departments of Pathology and Ophthalmology, Rush University Medical Center, 1 725 W. Harrison St., Suite 906, 

Chicago, Illinois 60612, USA 
6 Now with Lenticular Research Group, Naperville, Illinois 60540, USA 
*anastasio(a),seas. wustl. edu 



Abstract: The in vitro investigation of many optically opaque biological 
microstructures requires 3D analysis at high resolution over a large field of 
view. We demonstrate a new nondestructive volumetric imaging technique 
that eliminates the structural and computational limitations of conventional 
2D optical microscopy by combining x-ray phase-contrast tomography with 
critical point drying sample preparation. We experimentally demonstrate the 
enhancement of small features afforded by phase-contrast imaging and 
show the contrast improvement afforded by the drying of a hydrated 
specimen. We further demonstrate the biological application of this 
technique by imaging the microstructure of the accommodative apparatus in 
a primate eye using a benchtop phase -contrast tomography system. 

©2012 Optical Society of America 
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1. Introduction 

Although biological cells, tissues, and organs are three-dimensional (3D) structures, their in 
vitro physiology and pathophysiology are to this day primarily studied with two-dimensional 
(2D) imaging techniques (see, e.g., [1-3]). If the 3D nature of biological structures is not 
accounted for, key phenotypic and functional characteristics are often not recognized resulting 
in misinterpretations of the mechanisms responsible for both normal and disease processes. 

Analysis of microscopic 3D biological structures is not a simple task. In fact, the 
microscopic analysis of biological structures generally requires destructive manipulation of 
the 3D tissue structures that are to be studied. In order to examine tissue under a light or 
transmission electron microscope, for example, the tissue is generally embedded in wax or 
plastic and then sliced into sections. Although the sectioning process is well known to 
introduce structural artifacts [4,5], the process allows light or electron beams to pass through 
the tissue and thereby enables a 2D image of the sliced material to be captured. Only as the 
2D images of successive slices are layered together is the 3D nature of the tissue restored (see, 
e.g., [6]). Although computational techniques have facilitated and improved the accuracy of 
this process, a full accounting of 3D structures in intricate biological tissues remains a 
difficult, destructive, and time-consuming proposition. 

Overcoming these limitations is key to advancing our understanding of how biological 
form follows function and how malformation leads to malfunction. In this work, we describe a 
new nondestructive method that combines benchtop high-resolution phase-contrast x-ray 
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tomography with critical point drying to more fully account for microscopic 3D structures in 
opaque specimens. 

In x-ray phase-contrast tomography, variations in the spatial distribution of a sample's x- 
ray refractive index are employed to produce image contrast that augments or replaces 
conventional x-ray absorption contrast. Planar in-line phase-contrast imaging is a Gabor 
holography technique [7] that uses a highly-coherent input x-ray beam and a significant post- 
sample propagation distance to generate interference effects at tissue discontinuities. The 
resulting images contain a characteristic enhancement along sample edges and include 
intensity responses from small objects that can be practically invisible in traditional 
radiographic imaging. Gabor holography has been widely implemented with visible light 
sources [8,9]. In the x-ray regime, this technique was originally employed at synchrotron 
facilities [10-12] due to the brightness and high coherence of synchrotron x-ray sources, but 
the more recent use of microfocus x-ray tube sources, which produce fields with a high degree 
of spatial coherence, has enabled the creation of laboratory and clinical systems [13-15]. 

As in-line phase-contrast imaging derives contrast from refractive index variations, it is 
ideally suited for using air as a contrast agent. For example, studies have shown that the air- 
filled vasculature in liver specimens generates significant contrast by creating air-tissue 
refractive index discontinuities [16]. In an effort to study fine structures, such as the ciliary 
structures and zonules in the primate eye, we seek to exploit the contrast from air interfaces 
that may be introduced into these delicate microscopic structures. By marrying the critical 
point drying sample preparation technique, which is commonly used for electron microscopy, 
with x-ray phase-contrast tomography, the delicate structures of interest may be preserved 
while introducing contrast-rich air boundaries. The technique presented in this work allows 
for simultaneous large-scale and microscopic analysis without the destructive sectioning 
process that is required for microscopic evaluation. 

2. Materials and methods 

The benchtop in-line phase-contrast tomography system used in this work included a 
microfocus x-ray tube source (Thermo Kevex PXS10-65W; W anode, 7-100 um spot size, 65 
W maximum tube power, 45-130 kV peak voltage), a cooled high-resolution digital x-ray 
camera (Princeton Instruments Quad RO 4096; 16 urn pitch, 14-bit quantization), a three-axis 
sample positioning system (Thorlabs LTS150, XT95SP-1000), and a rotation stage for 
tomographic imaging (Thorlabs NR360S). The detector and sample positions were varied 
over the length of the table (2.4 m) to accommodate various specimen magnifications, system 
resolution requirements, and fields of view. In this study, the source-to-object distance R t 
varied from 14.3 cm to 68.4 cm and the object-to-detector distance R 2 varied from 25.2 cm to 
120.3 cm. The system spatial resolution was calculated as 



where s is the manufacturer-characterized spot size measured per the British Standard Method 
[17] and d « 33 um is the spatial detector blur estimated from a measured edge-spread 
function. Note that the tube current, which is directly related to the spot size and output flux 
from the microfocus tube source, was chosen to achieve the appropriate spot size and 
exposure time for a given experimental system geometry. The x-ray spectrum from the 
microfocus tube, when operated at 55 kV peak voltage, had an approximate minimum energy 
of 4 keV, a maximum energy of 55 keV, and an approximate mean energy of 16 keV. 

The imaging model for boundary-enhanced phase-contrast tomography is described in 
detail elsewhere [18]. Briefly, under a thin-object approximation, which is valid for tissue due 
to the high x-ray energy employed here, the measured intensity I(x, y) can be expressed in 
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(1) 



terms of the complex sample refractive index «(r) = l-£(r) + z'/?(r) , the system 
magnification M = (R l +R 2 )/ R t , and R 2 as 

I(x, y) = -^exp [-2* J &/?(r)] - ^ jexp [-2k\ dzfi(r)] J (2) 

where 7 0 is the unperturbed intensity with the object absent, k is the wave number, 

r = (x,y,z) is the position vector, z is the direction of propagation, and V^ >( is the 2D 

Laplacian operator acting on the detector-plane coordinates x and y. Note that the Laplacian 
operator acting on the right-hand term brings about edge-enhancement by acting on the real 
part of the projected complex refractive index, while conventional absorption contrast is 
brought about by the Beer-Lambert Law as described by the first term on the right-hand side 
of the equation. 

The function -ln(/(r) / I o ) is the projection data used for tomographic reconstruction of 

the 3D volume data via the Feldkamp-Davis-Kress (FDK) algorithm for circular cone -beam 
imaging [19]. The FDK algorithm is a direct reconstruction weighted convolution- 
backprojection formula that attempts to compensate for the finite cone-angle in contrast to 
slice-by-slice reconstruction techniques using parallel- or fan-beam data. The algorithm has 
been shown to be highly accurate for use in microfocus x-ray CT experiments. 

In this work, two specimens were selected for analysis. First, a formalin-fixed cranial 
defect [20] was used to demonstrate the phase-contrast enhancement brought about by the 
introduction of air into the sample via drying. This study was performed by removing the 
specimen from storage in a phosphate-buffered saline solution and drying via exposure to the 
room-temperature laboratory environment during imaging. Second, the anterior hemisphere of 
a previously-prepared critical point dried eye from an 11 -year-old baboon was used for 
tomographic analysis and to demonstrate the enhancement of fine structures. 

Critical point drying is commonly used to preserve small structures when preparing 
biological samples for scanning electron microscopy [21]. It avoids the creation of damaging 
surface tension forces associated with drying by bringing liquid in the sample to the gas phase 
without crossing the liquid-gas phase boundary. Briefly, the sample is generally prepared by 
flushing away water with acetone and then infusing it with high-pressure liquid carbon 
dioxide. The sample is subsequently heated until the pressure exceeds the critical point and 
the gaseous carbon dioxide is allowed to escape, leaving a dried sample. Great care must be 
taken when handling the dried specimen as fine structures are fragile. Detailed descriptions of 
critical point drying protocols are available elsewhere [22]. 

The critical point dried eye specimen was imaged using the benchtop in-line phase- 
contrast tomography system (Fig. 1). Tomographic data were acquired at 501 view angles 
(0.4-degree step) and reconstruction was performed using the FDK algorithm. Cornea and 
conjunctiva structures were removed from the reconstructed data by manual segmentation. 
Rendered 3D data were displayed using the Drishti volumetric visualization tool [23], which 
applies a two-dimensional transfer function to weight the opacity and color of each voxel. 

Changes in image appearance were quantitatively characterized by analysis of the image 
texture data. These data yield statistical information on the coarseness and smoothness of 
image content. Improved appearance of small structures, as is expected from phase -contrast 
imaging, results in increased coarseness of the image content, manifesting in increased sample 
variance and decreased uniformity. Image texture properties were quantified by analyzing the 
statistics of the image intensity histogram p(I) , where the overbar denotes a random variable 

and / is the intensity normalized to values of 0 to 255. The quantities of interest in this work 
are the sample variance 

255 

a 2 =^(I„-E) 2 p(I n ), (3) 

n=0 
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where E is the sample mean, and the uniformity 

255 _ 

71=0 

[24]. Low contrast, or smooth, regions are characterized by low sample variance and high 
uniformity. 

cooled 
detector 




Fig. 1. In-line phase-contrast tomography system. The planar x-ray intensity is shown as it 
propagates from the microfocus source to the sample (distance Ri) and from the sample to the 
detector (distance R 2 ). Intensity fringes develop with propagation and are measured by a digital 
detector. An example simulated phase-contrast image of a sphere is shown on the right. 

3. Results and discussion 

The introduction of air into the cranial defect sample is shown in Fig. 2. The specimen images 
were acquired with system parameters of source-to-object distance R ] =32.1 cm, object-to- 
detector distance R 2 = 65.8 cm, magnification M = 3.05, approximate spot size s » 13 urn, 

tube voltage V = 55 kV, and tube current J = 200 uA, yielding an effective pixel size of 4.9 
urn, and an estimated system resolution of r = 14 urn. As the specimen dried, the visibility of 
the of the fibrous poly-L-lactide scaffold microstructure (~40 urn fiber size) was improved 
dramatically. The change in contrast is due to the introduction of air interfaces into the 
specimen. As the refractive index difference between the specimen and surrounding air is 
much larger than the difference between the specimen and surrounding saline, the resulting 
refraction is much stronger and, thus, the phase contrast is improved. These results are 
qualitatively similar to those obtained in a previous study of air-contrast enhancement [14]. 

Changes to the histogram statistics are shown in Fig. 2, where each image was sub-divided 
into 25 regions for analysis and the calculated mean and variance (error bars) over the regions 
are shown. The increased sample variance (63% difference, from 1.3 x 10" 3 to 2.5 x 10" 3 ) and 
decreased uniformity (74% difference, from 7.9 x 10~ 2 to 3.6 x 10~ 2 ) with time indicates 
improved small feature contrast with drying. For comparison, an adjacent sample-absent 
image region containing only noise (not displayed) showed much smaller percentage 
differences in variance (18% difference, from 8.7 x 10" 3 to 1.0 x 10" 2 ) and uniformity (1% 
difference, from 0.121 to 0.122). The percentage differences are calculated as 
| of - a\ | /[(of + a\ ) / 2] and \U 1 -U 2 \ /[(£/, +U 2 )/2]. Finally, note that neither the phase- 
contrast technique nor the drying process limits the appearance of highly-absorbing calcified 
regions in these images. 

The enhancement of fine features due to phase-contrast imaging is shown in Fig. 3. The 
dried eye specimen was imaged ( V = 55 kV, J = 7 uA, s»5 urn) in two geometries, a low- 
propagation distance configuration (./?, = 14.3 cm, R 2 =25.2 cm) and a high-propagation 
distance configuration (R l =68.4 cm, R 2 - 120.3 cm). Note that the magnification is held 

constant (M - 2.76) in both cases so that the effective resolution of the detector is unchanged. 
It is clear from these data that the phase-contrast enhancement afforded by increased 
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Fig. 2. The contrast improvement induced by drying is shown here in projection images and 
plots of image histogram statistics. The cranial defect specimen, which was dried over a period 
of 58 minutes, shows fibrous polymer layer structures when dried (far right) that are relatively 
difficult to discern in the un-dried specimen (far left). Quantitative improvements in contrast 
are shown in the graphs of the time-dependent behavior of intensity histogram variance and 
uniformity. Note that the specimen occupies the entire frame in these images. The grayscale 
window is the same for all images. 



propagation distance improves the visualization of fine structures such as the layered structure 
of the crystalline lens. This result is expected, as increased propagation distance yields 
amplified edge -enhancement effects at discontinuities within the specimen. The resulting 
increased intensity variations in the high-propagation images were also expected to yield 
increased intensity variance and decreased intensity uniformity. 

Again, an analysis of image texture was performed to demonstrate the enhancement. The 
graph in Fig. 3 is a scatter plot showing the sample variance and uniformity over two regions 
in the low- and high-propagation images (indicated by boxes in Fig. 3(a) and Fig. 3(b)). Each 
region is divided into nine sub-regions, each denoted by a marker. For corresponding regions, 
the region-averaged variance is higher in the high-propagation distance image (1.7 x 10~ 2 ) 
than in the low-propagation distance image (2.9 x 10~ 3 , 141% difference) and the region- 
averaged uniformity is lower in the high-propagation distance image (1.7 x 10~ 2 ) than in the 
low-propagation distance image (4.2 x lfj -2 , 85% difference). This indicates that increased 
propagation distance increases the appearance of small features in the image data. For 
comparison, an adjacent sample-absent image region containing only noise showed much 
smaller percentage differences in variance (28% difference, from 7.2 x l(T 3 to 9.5 x 1(T 3 ) and 
uniformity (13% difference, from 0. 1 19 to 0. 104). 

In the computed tomography study ( V = 55 kV, J =400 uA, s « 30 urn, R t =45.2 cm, 

R 2 = 107.3 cm, M = 3.37), an estimated system resolution of r = 23 urn was achieved with 

effective pixel size 4.5 urn. The 3D tomographic data shown in Figs. 4 and 5 demonstrate the 
large field of view and high resolution of this technique. The 3D volume-rendered data 
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Fig. 3. The effect of propagation distance on phase contrast in projection images of a critical 
point dried baboon eye. In the low-propagation case (left), the fine structure of the crystalline 
lens and focusing structures are relatively indistinct. In the high-propagation case (right), in 
which phase-contrast is enhanced, these features are more easily visualized. Quantitative 
assessment of image statistics in two regions (boxes in the images) shows the effect of phase- 
enhancement on intensity histogram variance and uniformity. The grayscale window is the 
same for both images. 
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Fig. 4. Transverse slices of the reconstructed tomography data of a critical point dried baboon 
eye. The proximal plane slice in (a) clearly shows the iris (i), ciliary structures (ci), and choroid 
(ch). The distal plane slice in (b) additionally shows the crystalline lens (cl). Note that although 
there is little absorption contrast visible in the interior of many structures, boundaries display 
significant edge contrast (e.g., the iris). A threshold was applied to the data for display 
purposes. The grayscale window is the same for both images. 



(Fig. 5), which are shown in anterior, posterior, and lateral views, capture the entire 
accommodative apparatus of the eye, including the lens and ciliary structures. The fine 
focusing apparatus is best visualized in the zoomed 3D volume-rendered image shown in Fig. 
5(e), where the ciliary body, which appears as a set of fibrous structures surrounding the 
crystalline lens, is shown radiating outward from the lens body toward the posterior choroidal 
region. Also visible are the stratified lens structure and star-shaped surface sutures, which are 
well preserved during sample preparation and image reconstruction [25]. Transverse planar 
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slices (Fig. 4) show detailed images of the ciliary region including the stratified layers of the 
lens and the morphologic configuration of the lens, iris, and ciliary structures [26]. For 
reference, the iris thickness in adult humans has been measured to be approximately 0.5 mm 
and from high-magnification histology data, the thickness of the ciliary body is approximately 
0.2 mm [27,28]. 




Fig. 5. A 3D volumetric rendering of the focusing apparatus of a critical point dried baboon 
eye. Anterior (a), lateral (b), and posterior (c) images show the large field of view that can be 
captured relative to microscopy. The transverse slice data exposed in the anterior view (d) 
demonstrate the simultaneous visualization of interior and exterior sample structures. Finally, 
the zoomed-in posterior view (e) shows the highly detailed visualization of ciliary structures 
emanating from the crystalline lens. A threshold was applied to the volumetric data for display 
purposes. Scale bars are 1 mm in length. 

4. Conclusion 

X-ray phase-contrast tomography fills a significant void in the existing methods used for 3D 
analysis. Coupled with critical point drying, this imaging method holds great promise for a 
wide range of biological researches where improved methods for analysis of micro structure in 
optically opaque biological specimens are needed. Furthermore, the in-line phase-contrast 
tomography technique demonstrated here is widely applicable, as the laboratory system is 
contained on a benchtop, uses off-the-shelf components, and does not require access to an x- 
ray synchrotron source. 

In the eye research field, especially in the study of focusing structures, the methods 
presented here have a high potential impact. While electron microscopy techniques offer 
higher resolution, 3D imaging of an intact specimen is an untapped capability in the eye 
research community, which typically employs destructive 2D light or electron microscopy 
techniques to study morphology [1-3]. The method presented here is expected to augment 
these existing techniques to facilitate improved study of the 3D lens structure and 
accommodation. The latter topic is of particular interest as morphological changes in the fine 
structures making up the accommodative apparatus are being studied to gain an understanding 
of their effect on the focusing actions of the eye. These studies stand to be significantly 
improved by the simultaneous large-scale and microscopic specimen analysis described here. 
Moreover, by avoiding the specimen sectioning process required for light and electron 
microscopy, a significant source of structural defects has been eliminated. 
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In this work, we have shown that the combination of critical point drying and benchtop x- 
ray phase-contrast tomography yields highly detailed volumetric images of small delicate 
structures. This method holds significant promise as an analytic technique to augment 
microscopy and enable nondestructive evaluation of optically opaque specimens. 
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